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1.0  INTRODUCTION 

Stable,  aerodynamic  operation  of  the  compression  system  of  a  gas  turbine  engine  is 
essential  for  the  engine  to  perform  satisfactorily.  “Aerodynamically  stable  operation”  refers 
to  the  condition  existing  when  the  compressor  delivers  a  desired  quantity  of  airflow  on  a 
continuous  basis,  compressed  to  a  desired  pressure  ,  free  of  excessively  large  amplitude 
fluctuations  in  the  flow  properties  throughout  the  compressor.  The  aerodynamic  stability 
limit  (surge  line)— most  commonly  defined  as  a  locus  of  points  denoting  maximum  pressure 
ratio  for  stable  operation  as  a  function  of  corrected  airflow  rale — is  typically  determined 
from  experimental  results.  However,  only  limited  amounts  of  experimental  compressor 
steady-state,  transient,  and  stability  limit  data  can  be  obtained  during  test  programs  because 
of  test  hardware  and  economic  constraints.  A  validated  mathematical  model  based  on  the 
limited  experimental  data  is  needed  to  make  performance  and  stability  computations  at 
operating  conditions  with  inlet  distortion  patterns  and  with  subcomponent  hardware 
configurations  different  from  those  actually  tested.  The  use  of  such  a  model  would 
significantly  enhance  the  usefulness  of  data  from  any  given  compressor  test. 

The  objectives  of  the  work  described  herein  were  to  construct  a  validated  mathematical 
model  of  the  high-pressure  compressor  of  the  TF41-A-1  turbofan  engine  and  to  provide 
thereby  a  new  analysis  tool  for  assessing  the  influence  of  external  disturbances  and 
compressor  modifications  on  performance  and  stability.  The  objectives  were  accomplished 
through  development  of  a  one-dimensional,  steady-state  mathematical  model  for  stability 
assessment  with  undisturbed  flow  and  a  three-dimensional,  time-dependent  mathematical 
model  for  analysis  of  distorted  inflows  and  transient  and  dynamic  disturbances. 

For  both  models,  example  problems  are  presented  and  compared  to  experimental  results. 
The  problems  using  the  one-dimensional,  steady-state  model  consisted  of  determination  of 
the  steady-state  stability  limits  (surge  lines)  with  undisturbed  flow  for  three  distinct  inlet 
guide  vane  schedules.  Those  problems  using  the  three-dimensional,  time-dependent  model 
included  determination  of  the  stability  limit  (surge  line)  reduction  caused  by  pure  radial  inlet 
pressure  distortion,  pure  circumferential  inlet  pressure,  and  pure  circumferential  inlet 
temperature  distortion.  The  effects  of  rapid  upward  ramps  of  inlet  temperature  on 
compressor  stability  were  also  investigated.  The  models  compute  compressor  stability  limits 
with  reasonable  accuracy. 

An  altitude  evaluation  of  the  TF41-A-I  turbofan  engine  with  a  highly  instrumented, 
modified  (Block  76)  high-pressure  compressor  was  conducted  at  the  Arnold  Engineering 
Development  Center  during  the  period  September  15,  1977  through  January  26, 
1979.  Pressure  and  temperature  measurements  obtained  during  this  evaluation  provided 
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experimental  data  for  formulation  of  the  compressor  stage  characteristics  and  validation  of 
the  models. 

1 1 

The  technology  required  to  construct  compressor  models  applicable  to  predicting 
compressor  stability  margin  loss  caused  by  external  disturbance  in  Flow  conditions  was 
developed  in  Ref.  1.  This  demonstrated  the  validity  of  the  multidimensional  form  of  a 
mathematical  model  for  a  four-stage,  fixed -geometry  compressor  at  one  value  of  corrected 
rotor  speed.  In  this  program,  refinements  were  introduced  to  the  basic  (Ref.  1)  model.  These 
refinements  were  introduced  to  improve  the  generality  of  the  model,  to  accommodate  the 
variable  rotor  speed  and  the  variable-geometry  features  of  the  TF4I-A-1  high-pressure 
compressor,  to  increase  accuracy,  and  to  accommodate  the  inherently  unstable  numerical 
solutions  caused  by  error  buildup  in  the  model  calculations.  This  report  describes  the 
manner  in  which  the  models  were  constructed  and  the  results  obtained. 

2.0  TF41-A-1  TURBOFAN  ENGINE  ALTITUDE 
AEROMECHANICAL  TEST  AND  RESULTS 

Because  the  information  derived  from  the  TF41-A-1  turbofan  engine  altitude  test  is  such 
an  integral  part  of  the  effort  reported  herein,  a  description  of  the  engine  and 
instrumentation  and  pertinent  results  is  presented. 

2.1  ENGINE  DESCRIPTION 

The  TF41-A-1,  manufactured  by  the  Detroit  Diesel  Allison  (DDA)  Division  of  the 
General  Motors  Corporation,  is  a  two-spool,  low-bypass-ratio,  nonaugmented  lurbofan 
engine  (Fig.  1)  in  the  14,500-lbf-thrust  class.  The  engine  incorporates  a  three-stage  fan  and  a 
two-stage,  intermediate  compressor  driven  by  a  two-stage,  low-pressure  turbine  and  an 
eleven-stage,  high-pressure  compressor  driven  by  a  two-stage,  high-pressure  turbine.  The 
engine  incorporates  a  can-annular  combustor,  an  engine-controlled,  high-pressure 
compressor  interstage  (seventh-stage)  bleed  system,  and  an  annular  fan  duct.  A  schematic  of 
the  engine  installation  in  the  T4  test  cell  is  presented  in  Fig.  2. 

The  fan,  intermediate  compressor,  and  high-pressure  compressor  pressure  ratios  at  sea- 
level-static,  intermediate  power  are  2.48,  1 .41 ,  and  6.20,  respectively,  for  an  overall  pressure 
ratio  of  21.7  ar  a  rated  engine  airflow  of  260  lbm/sec,  a  high-pressure  rotor  speed,  NFL- 
of  12,840  rpm,  and  a  fan  rotor  speed,  NL,  of  8,940  rpm.  Airflow  is  split  at  the  fan  exit,  and. 
the  bypass  air  passes  through  the  annular  fan  duct  and  mixes  again  with  the  gas  generator 
airflow  aft  of  the  low-pressure  turbine  exit  at  the  tailpipe  inlet.  The  gas  generator  air  passes 
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through  the  intermediate  compressor  and  into  the  high-pressure  compressor  through  the 
variable-position  inlet  guide  vanes.  The  high-pressure  compressor  incorporates  a  bleed 
manifold  at  the  compressor  exit  in  the  outer  compressor  case  for  extraction  of  air  for  aircraft 
services.  The  fan  has  no  inlet  guide  vanes,  and  the  bulletnose  rotates;  there  are  no  anti-icing 
provisions. 

2.1.1  High-Pressure  Compressor  Modifications 

The  engine  (S/N  141683/3)  utilized  a  production  high-pressure  compressor  that  had  been 
modified  to  include:  (1)  full-chord,  clapperless,  first-stage  rotor  blades  with  increased 
thickness-to-chord  ratio  and  revised  taper,  (2)  “Eiffel  Tower”  stator  vanes  for  stages  4,  5, 
and  6  having  increased  thickness  near  the  root,  looser  dovetails,  and  glass  bead  peened 
surfaces,  and  (3)  modified  spacer  between  stages  9  and  10.  The  first  two  modifications  are 
referred  to  as  "Block  76”  design  modifications.  A  detailed  schematic  of  the  high-pressure 
compressor  is  presented  in  Fig.  3. 

2.1.2  High-Pressure  Compressor  Variable-Geometry  Control  System  1 

The  high-pressure  compressor  incorporates  a  hydromechanical  control  system  to 
position  the  inlet  guide  vanes  (IGV’s)  and  the  compressor  bleed  valve  (BVP)  by  a  common 
actuator  through  a  system  of  linkages. 

The  inlet  guide  vanes  are  scheduled  as  functions  of  engine  inlet  air  temperature  and  high- 
pressure  rotor  speed.  The  bleed  valve  is  fully  open  when  the  guide  vanes  are  at  35  deg  and 
fully  closed  when  the  guide  vanes  reach  24  deg. 

2.2  INSTRUMENTATION 

Instrumentation  was  provided  to  measure  dynamic  stresses  and  aerodynamic  pressures 
and  temperatures;  rotor  speeds;  fuel  system  pressures,  temperatures,  and  flow  rates;  lube  oil 
and  hydraulic  system  pressures  and  temperatures;  and  other  engine  and  test  cell  system 
parameters  required  for  proper  operation  of  the  engine. 

The  instrumentation  providing  data  for  this  study  was  located  in  the  high-pressure  (HP) 
compressor.  The  HP  compressor  was  heavily  instrumented  since  one  of  the  objectives  of  the 
altitude  lest  program  was  to  determine  the  detailed  aerodynamic  characteristics  of  the  HP 
compressor  throughout  the  TF41-A-1  engine  operating  envelope.  In  addition  to  the  array 
of  steady-state  instrumentation,  the  HP  compressor  was  also  instrumented  with  high- 
response,  wall-static  pressure  transducers.  A  diagram  showing  the  type,  number,  and 
location  of  the  HP  compressor  interstage  instrumentation  is  presented  in  Fig.  4. 
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2.3  PROCEDURE  FOR  HIGH-PRESSURE 

COMPRESSOR  STABILITY  ASSESSMENT 

Following  setting  of  the  desired  simulated  flight  conditions  the  high-pressure  rotor  speed 
was  adjusted  to  yield  a  selected  inlet  guide  vane  position,  IGV,  thereby  establishing  the 
baseline  point.  The  HP  compressor  was  then  slowly  back-pressured  by  opening  the  inbleed 
air  system  valve  until  surge  occurred;  steady-state  and  transient  data  were  obtained  along  the 
constant  IGV  line  from  baseline  to  the  surge  point.  The  power  lever  was  moved  during 
compressor  loading  to  maintain  the  selected  HP  IGV  angle  at  a  constant  value.  When  the 
HP  compressor  surged,  the  inbleed  air  system  valve  was  closed  to  relieve  the  compressor  and 
to  permit  recovery. 

The  surge  testing  of  interest  to  this  report  was  conducted  at  30,000  ft,  Mach  number  0.7 
in  the  range  of  operation  where  the  seventh-stage  bleed  valve  wras  closed  and,  at  sea-level- 
static,  where  open  seventh-stage  bleed  valve  performance  was  obtained.  To  evaluate  the 
effects  of  severely  off-schedule,  HP  compressor  IGV  geometry,  three  IGV  schedules  were 
selected,  representing  extremes  beyond  which  the  TF41-A-1  engine  is  not  expected  to 
operate.  These  IGV  schedules  are  defined  as: 

1.  “Nominal,”  M72-3  vane  schedule  (Ref.  2), 

2.  “Axial,”  open  approximately  10  deg  at  a  given  NH,  and 

3.  “Cambered,”  closed  approximately  10  deg  at  a  given  NH. 

These  schedules  are  shown  in  Fig.  5.  For  each  IGV  schedule,  a  HP  compressor  map  was 
obtained  by  inbleeding  high-pressure  air  into  the  exit  diffuser  by  the  previously  stated 
procedure. 

2.4  EXPERIMENTAL  RESULTS 

The  following  are  results  obtained  during  the  TF41-A-1  altitude  test  pertinent  to  this 
report.  The  normal  operating  and  surge  lines  for  the  HP  compressor  at  30,000  ft,  Mach 
number  0.7,  with  nominal,  axial,  and  cambered  IGV  schedules  are  presented  in  Fig.  6.  The 
compressor  surge  margin  was  26  percent  at  46.0  lbm/sec  corrected  airflow.  Misscheduling 
the  IGV  approximately  10  deg  axial  and  cambered  had  no  measurable  effect  on  the  Block  76 
compressor  surge  margin. 
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The  normal  operaiing  and  surge  lines  for  the  HP  compressor  at  sea-level-static 
conditions  with  nominal,  axial,  and  cambered  IGV  schedules  at  low  rotor  speeds  at  which 
the  seventh-stage  bleed  valve  is  open  are  presented  in  Fig.  7.  Compressor  exit  airflow 
corrected  to  compressor  inlet  conditions  was  used  as  the  abscissa  of  this  curve  because  no 
provision  was  made  to  calculate  seventh-stage  bleed  flow  rate,  and  high-pressure 
compressor  inlet  airflow  is  therefore  undefined  since  compressor  airflow  was  computed 
downstream  by  assuming  choked  flow  in  the  HP  turbine  vanes.  The  compressor  surge 
margin  with  the  bleed  valve  100  percent  open  was  41  percent  at  24.5  Ibm/sec  corrected 
airflow  with  the  nominal  IGV  schedule  and  at  27  Ibm/sec  with  the  cambered  IGV  schedule. 
Surge  could  not  be  induced  with  the  bleed  valve  100  percent  open  with  the  axial  IGV 
schedule. 

Interstage  dynamic  pressures  during  a  typical  surge  sequence  are  presented  in  Fig.  8. 
That  the  surge  originated  near  the  sixth  stage  is  indicated  by  the  pressure  traces  reflecting  the 
increase  in  dynamic  pressure  at  stages  I,  2,  and  4  and  by  the  corresponding  decrease  in 
dynamic  pressure  at  stages  8  and  10  following  the  initiation  of  surge. 

The  aerodynamic  pressures  and  temperatures  from  the  HP  compressor  vane  interstage 
instrumentation  (Fig.  4)  are  presented,  respectively,  in  Figs.  9  and  10  as  typical  examples  of 
the  total-pressure  and  total-temperature  profiles  existing  during  uniform  engine  inlet-flow 
conditions. 

Experimental  data  and  results  presented  are  compared  to  the  mathematical  model 
predictions  in  following  sections. 

3.0  DESCRIPTION  OF  ONE-DIMENSIONAL,  STEADY-STATE  MODEL 

A  one-dimensional,  steady-state  model  or  the  TF4I-A-1  HP  compressor  was  developed, 
based  on  compressor  geometry  provided  by  the  engine  manufacturer  and  on  stage 
characteristics  formulated  from  the  experiment  data.  The  steady-state  model  makes  use  of 
individual  stage  characteristics  of  pressure  and  temperature  rise  as  functions  of  mass  flow  to 
predict  ihe  quasi-steady  uniform  flow  performance  of  the  compressor. 

3.1  MODEL  CONCEPTS 

A  steady-state  compressor  model  is  simply  a  mathematical  representation  of  the 
compressor  performance.  Sufficient  data  are  necessary  to  determine  the  total-pressure  and 
total-temperature  ratios  for  the  compressor  as  functions  of  corrected  airflow  and  rotor 
speed.  Thus,  modeling  the  overall  steady-stale  performance  of  a  compressor  can  be  as 
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simple  as  reading  a  compressor  map.  However,  for  the  model  to  predict  stage  performance 
as  well  as  overall  performance,  it  is  necessary  to  analyze  the  components  or  stages  of  the 
compression  system  to  account  for  stage  interactions.  The  steps  for  development  of  the  one¬ 
dimensional,  steady-state  model  are  (1)  defining  control  volumes,  (2)  writing  governing 
equations,  and  (3)  giving  the  method  of  solution. 

3.1.1  Control  Volume  Analysis 

The  TF41-A-1  HP  compressor  system  (overall  control  volume)  was  divided  into  twenty 
elemental  control  volumes  consisting  of  four  inlet  duct,  eleven  compressor,  and  five  exit 
duct  elemental  control  volumes  (Fig.  11).  The  division  of  the  inlet  and  exit  ducting 
permitted  calculation  of  the  flow  properties  in  the  ducting  by  applying  conservation  of  mass, 
momentum,  and  energy  laws.  The  division  of  the  compressor  into  elemental  control 
volumes  representing  stage  lengths  in  the  axial  direction  was  made  because  the  pressure  and 
temperature  rise  between  stage  entry  and  exit  were  obtained  from  experimentally  determined 
stage  characteristics  as  discussed  in  Section  3.2.  For  the  one-dimensional,  steady-state 
analysis,  all  flow  properties  are  “bulked”  in  the  radial  and  circumferential  directions-,  and 
variations  with  the  axial  coordinate  only  are  considered. 

3.1.2  Governing  Equations 

An  axial-flow  compressor  imparls  kinetic  energy  to  the  flow  of  a  fluid  by  a  direct 
interaction  between  that  fluid  and  the  rotating  machine  elements.  The  fluid  is  then  diffused 
through  the  diverging  flow  passages  between  the  blades,  reducing  the  kinetic  energy  and 
increasing  the  static  pressure.  The  divergence  of  each  passage  is  determined  by  the  inlet  and 
exit  air  angles.  Since  the  blades  are  most  often  closely  spaced,  the  exit-air  angle  is 
approximately  equal  to  the  trailing-edge  angle  of  the  blades.  The  inlet-air  angle,  however,  is 
a  function  of  the  axial  velocity  entering  the  blades  and  their  rotational  speed.  Thus,  a  stage 
flow  coefficient  (<f>),  as  defined  in  Ref.  1, 


u 


where  U  =  axial  velocity  and  Uuh  =  wheel  speed,  was  used  to  represent  the  inlet-air  angle  at 
the  entry  of  each  compressor  elemental  control  volume. 

The  flow  coefficient,  0,  is  related  to  the  blade  angle  of  attack,  a,  by  the  stagger  angle,  X, 
which  is  a  constant  for  any  given  blade  row  through  the  relationship 
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<f>  ~  cat  (A  - «)  (2) 

Thus,  0  serves  the  same  purpose  in  the  compressor  stage  as  at  does  for  aircraft  wing 
aerodynamics.  It  decreases  with  increasing  or,  therefore,  small  0  means  large  a,  and  vice 
versa. 


The  temperature  rise  across  each  compressor  elemental  control  volume  was  represented 
by  the  stage  temperature  coefficient,  01 ,  expressed  as 


C  (Til— I) 

P 

uwt,2A 


(3) 


where  TR  is  the  stage  or  compressor  elemental  control  volume  temperature  ratio,  Cp  is 
specific  heat  at  constant  pressure,  and  T  is  the  stagnation  (total)  temperature. 


The  temperature  ratio  of  Eq.  (3)  for  an  ideal  compressor  stage  is  the  isentropic 
temperature  ratio  defined  from  stage  entry  and  exit  stagnation  pressure.  Using  this  ratio, 
another  stage  parameter,  the  stage  pressure  coefficient,  0P,  is  defined. 


(4) 


where  PR  is  the  stage  or  compressor  elemental  control  volume  pressure  ratio. 


cv  Cp  —  R/J  (5) 

and  where  y  is  the  ratio  of  specific  heats,  Cv  is  the  specific  heat  at  constant  volume,  and  R/I 
is  the  gas  constant/mechanical  equivalent  of  heat.  For  this  model,  the  specific  heat  at 
constant  pressure,  Cp,  was  computed  as  a  function  of  temperature  from  the  empirical 
equation: 


Cp  =  0.2318  h  0.104  x  10-*T  +  0.7166  x  10-BT2  (6) 

where  the  constants  are  based  on  the  specific  heats  of  the  constituents  of  air. 
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The  temperature  and  pressure  coefficients  (i/>T  and  ^p)  are  primarily  functions  of  the 
flow  coefficient,  (>,  and  the  flow  direction  (Ref.  1).  Thus,  a  set  of  curves 


=  f(*r 

^  =  f  (<«  . 


(7) 


is  referred  to  as  “stage  characteristics,”  and,  with  flow  direction  information,  fully  defines  a 
stage’s  performance  from  a  one-dimensional  standpoint. 


In  principle,  theoretical  or  experimental  lift  and  drag  coefficients  could  be  used,  along 
with  stage  geometry,  to  compute  the  stage  characteristics.  However,  corrections  for  blade 
row  interference  losses  and  secondary  flows  would  be  required.  The  method  used  in  this 
study  was  to  compute  directly  \fiT,  \pp,  and  <f>  by  Eqs.  (1),  (3),  and  (4)  from  the  experimental 
measurements  of  stage  total  temperature,  flow  rate,  and  total  pressure  at  the  stage  entry  and 
exit.  A  discussion  of  the  stage  characteristics  obtained  is  presented  in  Section  3.2. 


Friction  losses  in  the  inlet  and  exit  ducting,  internal  bleeds,  and  the  forces  acting  on  the 
fluid  in  the  compressor  must  satisfy  mass,  momentum,  and  energy  conservation  principles 
as  they  relate  to  any  elemental  control  volume  (Fig.  12). 

3.1. 2.1  Mass 


Application  of  the  mass  conservation  principle  to  the  elemental  control  volume  yields 


W  h—  dz  +  KG 

dr. 

”  V - 

mass  leaving 
control  volume 


Vi 


mass  entering 
control  volume 


(8) 


where  W  is  the  mass  flow  rate  and  z  is  the  axial  component. 
Equation  (8)  reduces  to 


W 


e  x  i  l 


w 

entrance 


-  WB 


(9) 


where  bleed  flow,  WB,  is  assumed  to  cross  the  control  volume  boundary  normal  to  the  axial 
direction. 
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3.1. 2. 2  Momentum 


The  summation  of  all  forces  in  a  given  direction  is  equal  to  the  change  of  momentum  in 
that  direction,  and  gives 


F  d*  +  PS  A  - 


6  (PS  A) 

A  + - 

dv 


axial  forces  acting  on  control  volume 


<9(WtI> 

dv 


momentum  leaving 
control  volume 


WU 

I  - - - M  1^ 

momentum  entering 
control  volume 


(10) 


where  F  is  the  Force  of  compressor  blading  and  cases  acting  on  fluid,  including  wall  pressure 
area  force;  PS  is  the  static  pressure;  and  A  is  the  area. 


Equation  (10)  may  be  reduced  to  give 

t?(IMP) 

$7. 


dA 

F  +  PS  — 


d. 


where  the  impulse  function,  IMP,  is  defined  as 


(II) 


IMP  =  WU  +  PS  A 


(12) 


The  resultant  axial  force,  F,  is  defined  as 

F  =  IMP  -  IMP 

' - v - ' 

impulse  leaving 
control  volume 

and  in  the  case  of  ducting  control  volumes,  assuming  adiabatic  flow  conditions,  reduces  to  a 
simple  friction  loss  expressed  as 


impulse  entering 
control  volume 


(13) 


F=-CdyPSAM2  (14) 

where  Cd  is  the  duct  pressure  loss  coefficient  and  M  is  the  Mach  number.  The  value  of  Cj 
was  adjusted  to  give  a  reasonable  stagnation  pressure  loss  in  the  ducting. 
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3. 1.2.3  Energy 


Energy  conservation  gives 


dll 

II  h - dz  = 

H 

+  WS  d  z 

t  J 

Qdz 

L  .  J 

- v - — ' 

V 

enthalpy  leaving 

enthalpy  entering 

shaft  work 

heat  added  to  fluid 

control  volume 

control  volume 

done  on  fluid 
in  control  volume 

in  control  volume 

(15) 

where  H  is  the  total  enthalpy  flux,  WS  is  the  stage  staft  work  added  to  fluid  in  control 
volume,  and  Q  is  the  rate  of  heat  addition  to  control  volume.  Equation  (15)  may  be  reduced 
to 


^  =  WS  +  Q 

dz  (16) 

where 

H  =  wcpT  (17) 

The  shaft  work,  WS,  is  defined  as 

WS  =  WC  T  (TR  - 1)  (18) 

and  is  zero  in  the  ducting  control  volumes. 

3. 1.2.4  State  and  Additional  Equations 

Additional  equations  required  include  the  perfect-gas  equation  of  state, 

PS  =  p  R  TS,  "  (19) 

where  p  is  the  density  and  R  is  the  gas  constant;  stagnation  (total)  state  temperature 
equation, 

li2 

T  =  TS  +  —  (20) 

2  Cp 
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and  stagnation  (total)  state  pressure  equation. 


y 


Also  required  is  the  Mach  number. 


u 

M  =  - 


where  the  acoustic  velocity,  a,  is 


(21) 


(22) 


(23) 


3.1.3  Method  of  Solution 

The  one-dimensional,  steady-state  model  uses  a  “stacking”  procedure  to  determine  flow 
properties  throughout  the  compressor  and  associated  ducting.  The  one-dimensional  model 
of  the  TF41-A-1  HP  compression  system  uses  segmented  inlet  and  exit  ducting  as  indicated 
in  Fig.  11.  The  solution  is  begun  with  known  air  properties  at  the  compressor  inlet  (total 
pressure,  total  temperature,  and  airflow),  ducting  geometry,  and  a  coefficient  of  friction. 
Flow  properties  throughout  the  inlet  ducting  are  calculated  using  the  one-dimensional, 
steady-state  conservation  laws  (see  Section  3.1.2).  Since  the  exit  of  one  duct  segment  is  the 
entrance  to  the  next  segment,  the  solution  is  marched  upstream  from  the  compressor  inlet  to 
the  duct  inlet  by  an  iterative  procedure.  For  the  compressor  segments,  with  compressor  inlet 
flow  properties  known,  a  stage-by-stage  calculation  downstream  through  the  compressor 
using  the  stage  characteristics  formulated  from  experimental  data  gives  individual  stage 
pressure  and  temperature  ratios  so  that  overall  compressor  performance  can  be  calculated. 
Compressor  aft-ducting  flow  properties  are  then  calculated  using  the  one-dimensional, 
steady-state  conservation  laws. 

3.2  TF41-A-1  HIGH-PRESSURE  COMPRESSOR  STEADY-STATE 
STAGE  CHARACTERISTICS 

Representative  stage  characteristics  (temperature  and  pressure  coefficients  as  functions 
of  flow  coefficient)  are  presented  schematically  in  Fig.  13  for  a  fixed-geometry  compressor 
operating  at  a  single  value  of  compressor  corrected  rotor  speed. 
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Several  problems  had  to  be  resolved  because  of  the  variable  geometry  of  the  TF41-A-1 
HP  compressor: 

1.  To  allow  investigation  of  off-nominal  schedules,  the  model  had  to  account  for 
variable  IGV  effects  without  relating  IGV  position  to  rotor  speed. 

2.  To  determine  the  stage  flow  rates  during  bleed-valve  open  operation  the  seventh- 
stage  bleed  flow  rate  had  to  be  estimated. 

3.  Interpolation  of  stage  characteristics  between  lines  of  the  chosen  dependent 
variables  had  to  reflect  the  experimentally  determined  results. 

Assumptions  were  made,  and  procedures  were  established  and  incorporated  into  the 
TF41-A-1  HP  compressor  models  to  meet  these  requirements.  The  steady-state  stage 
characteristics  are  utilized  both  in  the  one-dimensional,  steady-state  model  and  in  the  three- 
dimensional,  time-dependent  model;  therefore,  the  following  discussion  will  apply  to  both 
models. 

The  variation  of  the  steady-state  stage  characteristics  with  rotor  speed  was  assumed  to  be 
second  order,  with  primary  effect  attributed  to  the  IGV  position.  This  relates  the  stage 
characteristics  to  the  single  dependent  variable,  the  IGV,  and  greatly  simplifies  their 
formulation. 

An  empirical  relationship  between  bleed  valve  position,  BVP,  and  bleed  flow  rate,  WB, 
was  developed,  based  on  the  reactions  of  various  static-to-total-pressure  ratios  determined 
from  the  interstage  instrumentation  (Fig.  4).  This  estimation  of  the  seventh-stage  bleed  flow 
rate  permitted  calculation  of  the  stage  characteristics  from  experimental  data  during  open 
bleed-valve  operation  (see  Section  2.1.2). 

The  stage  characteristic  parameters  ^T,  i/,p,  and  4>  were  computed  from  measurements  of 
rotor  speed,  flow  rate,  total  pressure,  and  total  temperature  for  each  stage  entry  by  use  of 
Eqs.  (3),  (4),  and  (1)  respectively  for  all  data  obtained  during  the  TF41-A-1  test.  The 
stage  characteristics  were  then  generalized  along  lines  of  constant  IGV  angle  over  the  full 
range  of  IGV  operation  (from  -5.5  to  36  deg).  The  TF41-A-1  HP  compressor  stage 
characteristics  developed  from  this  analysis  are  presented  in  Fig.  14. 

The  stage  characteristics  were  incorporated  by  means  of  a  series  of  second-order  curve 
fits.  Frequently,  a  prediction  of  compressor  performance  is  required  at  some  value  of  IGV 
that  is  not  equal  to  any  of  the  values  used  to  formulate  the  stage  characteristics.  Hence, 
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interpolation  of  the  required  parameter  (f  or  ^p)  from  the  values  of  stage  characteristics 
between  higher  and  lower  guide  vane  positions  is  required.  The  interpolation  procedure  used 
in  the  models  is  based  on  a  fixed  value  of  <t>.  This  procedure  was  adopted  because  the 
variation  of  compressor  mass  flow  with  rotor  speed  is  effectively  represented  by  4>.  Hence, 
the  same  </>  can  be  used  to  determine  or  at  the  higher  and  lower  guide  vane  positions. 
The  i p1  or  ^p  can  then  be  interpolated  linearly  to  define  the  value  at  the  desired  IGV  for  a 
given  value  of  d>. 

3.3  MODEL  VERIFICATION 

Figure  15  presents  steady-state  model  computations  of  compressor  performance  at 
30,000  ft  and  Mach  number  0.7,  with  nominal,  cambered,  and  axial  inlet  guide  vane 
schedules  over  the  corrected  rotor-speed  range  from  90  to  100  percent  and  compares  them  to 
experimental  results.  The  steady-state  model  computes  compressor  surge  pressure 
ratios  within  2.4  percent  of  the  experimental  stability  limits  for  the  nominal,  cambered,  and 
axial  guide  vane  schedules.  For  these  results  (Fig.  15)  the  model  computations  represent 
constant  corrected  rotor-speed  lines  at  a  constant  IGV,  whereas  the  experimental  data  were 
obtained  with  nearly  constant  IGV  (  +  0.4  deg)  while  corrected  HP  rotor  speed,  NH/^ 
varied  as  much  as  3.6  percent  .The  procedure  used  to  obtain  the  experimental  data  (see 
Section  2.3),  by  maintaining  constant  IGV  during  the  HP  compressor  loading,  resulted  in 
near  constant  NH  (±1.0  percent).  However,  the  increase  in  energy  to  the  engine  turbines 
resulted  in  increased  fan  rotor  speed,  NL.  Thus,  the  inlet  temperature,  T24,  delivered  to  the 
HP  compressor  increased  from  the  baseline  point  to  the  stall  point  with  corresponding 
decreases  in  NHA/0-m.  Therefore,  the  model  and  experimental  data  represent  different  paths 
from  baseline  to  stall  point.  This  is  illustrated  in  Fig.  16  in  which  the  steady-state  model 
predictions  are  presented  at  the  corresponding  corrected  rotor-speeds  of  the  experimental 
data  points  for  the  -0.3-deg  IGV  line  on  the  axial  IGV  schedule  (Fig.  15c). 

Steady-state  model  computations  of  overall  compressor  performance  with  seventh-stage 
bleed  valve  open  at  sea-level-static  conditions  with  nominal,  cambered,  and  axial  inlet  guide 
vane  schedules,  over  the  corrected  rotor-speed  range  from  85  to  92  percent  are  presented  and 
compared  to  experimental  results  in  Fig.  17.  The  steady-state  model  computes  compressor 
surge  pressure  ratios  within  4.0  percent  of  the  experimental  stability  limits  for  the 
nominal,  cambered,  and  axial  guide  vane  schedules  with  seventh-stage  bleed  valve  open. 

The  large  quantity  of  steady-stale  temperature  and  pressure  measurements  made  on  the 
TF41-A-1  HP  compressor  during  testing  allows  detailed  model-to-experiment 
comparisons  on  a  stage-by-stage  basis.  Figures  18,  19,  and  20  show  model  computed  and 
experimentally  determined  stage  flow  coefficients,  pressure  ratios,  and  temperature  ratios. 
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respectively,  as  functions  of  corrected  rotor  speed  near  the  compressor  stability  limit.  Good 
agreement  was  obtained  (within  0.5  percent). 

It  is  evident  from  Fig.  18  that  theTF41-A-l  HP  compressor  operates  with  forward  stages 
(1  through  3)  below  the  stall  value  of  stage  flow  coefficient  at  corrected  rotor  speeds  below 
approximately  91  percent.  This  operation  reinforces  previous  theories  and  experiment 
observations  (Refs.  1 , 3,  4,  and  5).  At  low  corrected  rotor  speeds,  the  forward  stages  cannot 
deliver  their  intended  pressure  ratio.  The  flow  density  provided  to  the  aft  stages  is  therefore 
lower  than  that  for  which  the  compressor  axial  area  distribution  was  designed.  Thus,  to 
maintain  steady  mass  flow  continuity,  the  flow  coefficients  (axial  velocity)  in  the  aft  stages 
must  increase  with  decreasing  corrected  rotor  speed.  Indeed,  at  sufficiently  low  corrected 
rotor  speeds,  the  aft  stages  are  even  choked.  At  low  corrected  rotor  speeds,  as  the  load  on 
the  compressor  (the  demanded  pressure  ratio)  is  increased  by  back-pressuring,  the 
compressor  flow  reduces,  and  the  following  occur:  Initially  the  aft  stages  may  increase  their 
pressure  ratio  in  response  to  the  flow  reduction  at  a  greater  rate  than  the  rate  of  decrease  in 
pressure  ratio  occurring  in  the  forward  stages;  thus,  the  increased  load  may  be  satisfied. 
However,  as  the  loading  continues,  the  forward  stages  become  more  and  more  stalled,  and  a 
point  is  eventually  reached  at  which  the  aft  stages  can  no  longer  compensate.  At  that  point, 
overall  compressor  instability  occurs. 

It  is  necessary  to  distinguish  stage  stalling  from  overall  compressor  stalling  or 
instability.  Instability  is  always  preceded  by  stage  stall  but  is  not  always  simultaneous  with 
stage  stall. 


4.0  DESCRIPTION  OF  THREE-DIMENSIONAL, 
TIME-DEPENDENT  MODEL 

The  three-dimensional,  lime-dependent  model  detailed  in  Ref.  1  was  converted  to  the 
TF41-A-1  high-pressure  compressor  system  by  the  incorporation  of  TF41-A-1  compressor 
geometry  and  stage  characteristics  formulated  during  the  steady-stale  model  investigation 
discussed  in  Section  3.  The  lime-dependent  model  utilizes  a  finite  control  volume,  integral 
approach  to  solve  the  full  nonlinear  form  of  the  three-dimensional  conservation  laws  of 
mass,  momentum,  and  energy  by  means  of  a  finite  difference  technique.  Thus,  the  model  is 
capable  of  computing  the  influences  of  transient  and  dynamic  disturbances  with  either 

uniform  or  distorted  inflow. 

( 

4.1  MODEL  CONCEPTS 

A  detailed  description  of  the  concepts  involved  in  the  three-dimensional,  time-dependent 
model  development  is  presented  in  Ref.  1  and  summarized  in  this  section. 
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The  compressor  was  divided  radially,  circumferentially,  and  axially  on  a  stage-by-stage 
basis  into  control  volumes  (Figs.  21  and  22).  The  governing  equations  were  developed  from 
the  three-dimensional,  time-dependent  mass,  momentum,  and  energy  equations  for  a  finite 
control  volume.  The  force  and  shaft  work  applied  to  the  fluid  are  determined  from  empirical 
steady-state  stage  characteristics  modified  for  unsteady  cascade  airfoil  effects.  For 
simplification  empirical  crossflow  relationships  were  used  to  replace  the  radial  and 
circumferential  momentum  equations.  The  radial  work  distribution  variation  effect 
necessary  for  treatment  of  radial  distortion  was  approximated  using  average  stage 
characteristics  modified  as  a  function  of  radius  by  empirical  total  pressure  and  total 
temperature  profile  correction.  Crossflow  and  unsteady  cascade  effects  for  circumferential 
distortion  were  also  built  into  the  model. 


4.1.1  Method  of  Solution 


The  three-dimensional,  time-dependent  compressor  model  is  formulated  as  an  initial 
condition  problem.  The  values  of  the  dependent  variables  (mass,  momentum,  and  energy) 
are  specified  for  every  control  volume  location  (axial,  radial,  and  circumferential 
coordinate)  to  start  the  solution.  For  simplicity,  the  solution  is  always  started  from  a  steady, 
uniform  flow  condition,  thus  allowing  determination  of  the  initial  conditions  by  use  of  the 
one-dimensional,  steady-state  model.  With  dependent  variables,  p,  pU,  and  X  (representing 
mass,  momentum,  and  energy)  specified  for  every  control  volume  location,  ijk,  for  time 
equal  zero,  the  time  derivatives  are  computed  from  the  following  equations: 


4.1. 1.1  Mass 
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4. 1.1. 2  Axial  Momentum 
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The  force,  F,  is  obtained  from  stage  characteristic  information  modified  for  unsteady 
cascade  airfoil  effects.  The  impulse  terms  are 


IMP  =  PS  AZ  +  WZ  U 


(26) 


where  AZ  is  the  control  volume  area  normal  to  axial  direction  and  WZ  is  the  mass  flux 
across  radially-facing  control  volume  boundary. 
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4. 1.1.3  Energy 
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where  WS  is  the  stage  shaft  work  added  to  fluid  in  control  volume,  Q  is  the  rate  of  heat 
addition  to  control  volume,  HZ  is  the  total  enthalpy  flux  transported  across  control  volume 
axial-facing  boundary,  and  where 


X 


=  P 


(28) 


The  radial  and  circumferential  impulse  and  velocity  terms  neglected  in  Eqs.  (25)  and  (28)  are 
second-order  terms  determined  respectively  by  nondimensionalization  and  order-of- 
magnitude  analysis  of  the  equations;  therefore,  their  influence  on  the  solution  would  have 
been  small.  Also,  including  these  terms  caused  the  numerical  solutions  to  be  generally  less 
stable  and  to  oscillate  with  small  amplitude  when  a  steady-state  condition  was  reached.  Once 
these  calculations  of  the  time  derivatives  of  the  dependent  variables,  p,  pU,  and  X,  are 
carried  out  for  all  the  volumes,  the  solution  can  be  advanced  in  time  through  use  of  a  fourth- 
order  Runge-Kutta  numerical  integration  method  (Ref.  6)  that  provides  values  of  the 
variables  at  the  next  time  step.  The  sequence  is  then  repeated,  with  boundary  conditions 
changing  in  accordance  with  the  specified  event  simulation  for  as  many  lime  steps  as  are 
required.  A  schematic  of  the  computational  process  is  presented  in  Fig.  23. 


The  system  of  equations  constituting  the  model  is  for  time-dependent,  three-dimensional 
subsonic  flow  and  is  hyperbolic  in  nature,  with  specified  boundary  conditions  and  initial 
conditions.  Thus,  each  integration  step  to  compute  the  dependent  variables  at  the  new  time 
(t  +  At)  must  be  made  within  the  region  of  influence  of  the  known  values  at  time,  t:  The 
Courant-Friedrichs-Lewy  maximum  time  step  criterion  (Ref.  7) 
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must  be  observed.  For  AZmin  =  0.111  ft,  a  =  1,625  ft/sec  (values  at  TS  =  1,100°R), 
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A  time  step  of  1.0  x  lO-5  sec  was  used  for  the  TF41-A-1  HP  compressor  model. 

4.2  THREE-DIMENSIONAL,  TIME-DEPENDENT  MODEL  REFINEMENTS 

Before  this  study,  the  three-dimensional,  time-dependent  model  concepts  (Ref.  1) 
summarized  in  the  previous  section  had  only  been  applied  to  a  four-stage,  fixed-geometry 
compressor  operating  at  a  single  rotor  speed.  During  this  investigation,  a  number  of 
refinements  were  introduced  to  the  basic  model  to  accommodate  variable  rotor  speed  and 
variable-geometry  features  of  the  TF41-A-1  HP  compressor.  In  addition,  refinements  were 
incorporated  to  improve  accuracy,  and  to  handle  numerical  solutions  that  were  inherently 
unstable.  (This  instability  was  caused  by  error  buildup  in  the  model  calculations.)  These 
refinements  are  discussed  in  the  following  paragraphs. 

4.2.1  Variable  Rotor  Speed  and  Variable  Geometry 

The  variable-rotor-speed  and  variable  geometry  requirements  were  met  by  incorporating 
the  IGV-dependent,  steady-slate  stage  characteristics  formulated  during  the  one- 
dimensional,  steady-state  model  investigation  (Section  3.2).  With  both  NH  and  1GV 
specified  as  inputs  (boundary  conditions),  stage  shaft  work  and  forces  could  be  determined 
for  all  1GV  schedules. 

4.2.2  Real-Gas  Properties 

As  was  done  in  the  one-dimensional,  steady-state  model  development,  real-gas 
properties  were  approximated  by  incorporation  of  an  emperical  relationship  for  the  specific 
heat  at  constant  pressure,  Cp,  as  a  function  or  temperature  and  gas  composition. 

4.2.3  Modified  Fourth-Order  Rungc-Kulla  Numerical  Integration  Method 

Conversion  of  the  three-dimensional,  time-dependent  model  to  the  TF4I-A-1  HP 
compressor  form  with  incorporation  of  the  previously  discussed  refinements  resulted  in  a 
mathematical  model  that  was  numerically  unstable  at  all  conditions.  All  attempts  to  start  the 
model  resulted  in  unsteady  conditions  that  produced  a  model  computation  abort  similar  in 
nature  to  a  compressor  surge.  An  analysis  of  the  instability  revealed  that  the  fourth-order 
Runge-Kutta  numerical  integration  method  was  not  convergent  when  applied  to  the 
hyperbolic  equations  constituting  the  TF41-A-1  HP  compressor  model.  Further  analysis 
revealed  that  the  divergence  was  of  a  series  form,  alternating  around  the  desired  solution. 
Therefore,  an  average  of  the  oscillations  over  at  least  one  cycle  approximates  the  desired 
solution. 
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A  technique  was  developed  and  incorporated  to  average  the  oscillations  of  the  dependent 
variables  (mass,  momentum,  and  energy)  over  a  predetermined  time  interval  and  restarting 
the  solution  at  the  average  value  of  the  dependent  variables.  The  averaging  technique  was 
investigated  for  convergence  by  running  various  time  intervals  to  determine  the  minimum 
interval  at  which  convergence  could  be  maintained.  The  minimum  time  interval  for 
satisfactory  convergence  was  determined  to  be  approximately  0.2  x  10'  sec.  This  time 
interval  was  also  effective  at  the  time  of  actual  compressor  surge  since  the  instabilities  that 
occur  during  the  surge  event  fully  develop  and  produce  model  computation  abort  within  0.2 
x  10-3  sec.  An  illustration  of  the  effect  of  the  incremental  averaging  technique  on  numerical 
instability  is  presented  in  Fig.  24. 

4.3  MODEL  VERIFICATION 

The  TF41-A-1  HP  compressor  model  was  loaded  to  the  stability  limit  by  increasing  the 
exit  static  pressure  until  flow  breakdown  was  indicated.  Inlet  stagnation  pressure  and 
temperature  were  constant.  Corrected  rotor  speed  was  92.6  percent  of  design,  and  IGV  angle 
was  12  deg  (nominal  IGV  schedule).  Figure  25  illustrates  the  loading  and  the  resulting  stall 
process.  Exit  static  pressure  was  ramped  at  a  rapid  rate  to  minimize  the  computer  time 
necessary  for  the  computations.  The  rate  was  slow  enough,  however,  that  transient  effects 
were  very  small.  Airflow  at  any  stage  never  differed  from  compressor  inlet  airflow  by  more 
than  0.6  percent  during  loading.  As  the  load  on  the  compressor  increased,  compressor 
pressure  ratio  increased,  and  airflow  decreased  until  flow  breakdown  was  indicated. 

Figure  26  presents  the  model  computed  stage  inlet  total  pressures  on  an  expanded  time 
scale  covering  the  flow  breakdown.  Initiation  of  surge  is  indicated  at  t  =  76.78  msec  by 
observation  of  the  sixth-  and  seventh-stage  inlet  total  pressures.  The  inlet  total  pressure  at 
stage  seven  decreases  because  stage  six  stalls,  thereby  limiting  pumping  capacity.  The  inlet 
total  pressure  of  stage  six  increases  in  response  to  the  sudden  blockage  to  the  flow  produced 
by  the  stalling  in  stage  six.  The  indication  that  compressor  surge  originated  in  the  sixth  stage 
agrees  with  the  experimental  results  reported  in  Section  2.4 

.  The  model  solution  thus  indicated  the  following  sequence  leading  to  compressor 
instability:  The  compressor  is  loaded  and  mass  flow  decreases  until  a  stage  somewhere  in  the 
compressor  reaches  its  stall  point  (its  maximum  pressure  coefficient  point).  The  ability  of 
that  stage  to  continue  to  pump  flow  (i.e.,  to  help  support  the  pressure  gradient  over  the 
length  of  the  compressor  imposed  by  the  back  pressure  load)  begins  to  diminish.  As  the  flow 
reduction  continues,  the  ability  of  the  stage  to  pump  diminishes  even  further  because  of  the 
increased  load.  Eventually,  a  point  is  reached  where  the  flow  into  the  stage  becomes  fully 
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stalled;  thus,  little  force  on  the  fluid  is  provided  to  pump  the  flow,  and  the  flow  path 
through  the  compressor  becomes  blocked  because  of  the  highly  separated  condition  of  the 
stalled  flow  in  the  channels. 

\ 

In  terms  of  the  solution  of  the  model  equations,  as  the  loading  is  increased,  a  point  is 
reached  where  no  stable  solution  to  the.  governing  equations  exists  that  will  satisfy  the 
imposed  load  (boundary  conditions);  thus,  the  solution,  like  the  flow  through  the 
compressor,  becomes  unstable.  However,  the  behavior  of  the  model  quickly  becomes 
nonrepresentative  of  the  ’  compressor  shortly  after  flow  breakdown  because  the  stage 
characteristics  are  not  valid  for  deeply  stalled,  dynamic  situations. 

4.3.1  Pure  Radial  and  Circumferential  Distortion  Patterns 

Three  types  of  distortion  patterns  were  computed  to  see  the  influence  of  pattern  severity 
and  shape  variations  on  the  TF41-A-1  HP  compressor.  The  patterns  considered  were  pure 
radial  pressure  distortion,  pure  circumferential  pressure  distortion,  and  pure  circumferential 
temperature  distortion.  No  specific  experimental  data  for  the  TF41-A-1  HP  compressor 
were  available  for  these  cases.  Therefore,  only  comparisons  to  experimental  results  from 
other  compressors  could  be  made. 

4.3.1. 1  Radial  Pressure  Distortion 

Figure  27a  illustrates  the  results  of  loading  (by  increasing  exit  static  pressure)  the  TF4I- 
A-l  HP  compressor  model  to  the  stability  limit  with  tip  radial  distortion  (low  pressure  in  the 
tip  region).  Pattern  severities  (Pmax  -  Pmii/Pavg).  of  5,  10,  and  20  percent  were  considered. 
Figure  27b  shows  the  results  for  the  hub  radial  case  and  indicates  that  the  TF41-A-1  HP 
compressor  would  be  more  sensitive  to  tip  radial  distortion  than  to  hub  radial  distortion. 
Figure  28  compares  of  the  TF41-A-1  HP  compressor  model  results  with  experimental  data 
from  a  J85-13  engine  test  (Ref.  8).  Simitar  results  were  obtained,  except  that  the  J85-13 
engine  was  even  less  tolerant  of  tip  radial  distortion  than  the  TF41-A-1  HP  compressor  was 
computed  to  be. 

4.3. 1.2  Circumferential  Pressure  Distortion 

The  effect  of  increasing  the  severity  of  a  180-deg-arc  circumferential  pressure  distortion 
pattern  is  presented  in  Fig.  29.  Compressor  surge  limit  in  terms  of  pressure  ratio  decreased 
3.1, 4.8,  and  7.9  percent,  respectively,  from  the  uniform  flow  (no  distortion)  surge  limit  at 
pattern  severities  of  5,  10,  and  20  percent.  Figure  30  shows  the  TF41-A-1  HP  compressor 
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model  results  compared  to  a  range  of  experimental  data  from  other  compressors  (Refs.  8 
and  9).  The  trend  and  magnitude  predicted  by  the  model  agree  well  with  the  experimental 
data. 

4.3. 1.3  Circumferential  Temperature  Distortion 

Figure  31  shows  the  effect  of  varying  the  severity  of  180-deg-arc  total  temperature 
distortion  on  compressor  stability.  Temperature  distortion  destabilizes  the  compressor  in 
much  the  same  manner  as  does  pressure  distortion.  The  higher  temperature  region  operates 
at  reduced  airflow,  thus,  at  reduced  flow  coefficients  compared  to  the  low  temperature 
region  of  the  compressor.  As  the  load  on  the  compressor  is  increased,  the  high  temperature 
region  stalls  first,  eventually  causing  the  entire  compressor  to  become  unstable. 

For  the  example  presented  in  Fig.  31,  corrected  rotor  speed  was  held  constant  during  the 
loading  (exit  static  pressure  increase)  by  maintaining  constant  average  temperature  over  the 
face  of  the  compressor. 

Figure  32  shows  a  comparison  of  the  TF41-A-I  HP  compressor  model  computed  stability 
limit  pressure  ratio  reduction  caused  by  circumferential  temperature  distortion  to  test  results 
with  other  compressors  (Refs.  8  and  10).  The  trend  and  magnitude  predicted  by  the  model 
agree  well  with  the, experimental  data. 

i 

4.3.2  Rapid  Inlet  Temperature  Ramp  to  Stall 

The  TF41-A-1  HP  compressor  model  was  subjected  to  a  series  of  uniform  inlet 
temperature  ramp  rates  of  2,000,  5,000  and  10,000° R/sec  to  simulate  ingestion  of  hot  gases 
from  gun  or  rocket  fire.  The  effects  of  inlet  temperature  ramps  were  investigated  for  three 
IGV  schedules  of  the  TF41-A-I  HP  compressor.  Guide  vane  schedules  investigated  were 
nominal,  cambered  ( +  10  deg  from  nominal),  and  axial  (-10  deg  from  nominal).  Mechanical 
rotor  speed  and  inlet  total  pressure  were  held  constant  while  exit  static  pressure  was 
maintained  at  the  level  corresponding  to  experimentally  determined  operating  line  condition 
for  the  respective  guide  vane  angle  investigated.  The  initial  corrected  rotor  speed  was  92.6 
percent.  Figure  33  shows  the  effect  of  5,000° R/sec  uniform  inlet  temperature  ramps  with  the 
nominal,  cambered  (+  10  deg  from  nominal)';  and  axial  (-10  deg  from  nominal),  mechanical 
no  significant  changes  in  transient  surge  limit  with  IGV  schedule;  this  is  expected  since  both 
steady-state  model  and  experimental  results  indicated  that  misscheduling  the  inlet 
guide  vane  had  no  measurable  effect  on  the  TF41-A-1  compressor  steady-state  surge  margin. 
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Analysis  of  the  model  computations  (Fig.  34)  indicate  that  the  forward  stages  arc  rapidly 
driven  into  stall  by  the  reduction  in  airflow  accompanying  the  rapid  temperature  rise  at  the 
compressor  inlet.  Initiation  of  surge  and  associated  flow  breakdown  occured  in  the  second 
stage  and  propagated  back  through  the  compressor. 

Figure  35  shows  the  TF4I-A-1  HP  compressor  model  computed  influence  of  uniform 
inlet  temperature  ramp  rates  on  compressor  stability  with  nominal,  cambered  and  axial  IGV 
schedules.  The  model  results  from  Fig.  36  indicate  that  the  nominal  schedule  produces 
slightly  less  reduction  in  airflow  for  a  given  temperature  rise  at  the  inlet  than  the  cambered 
or  axial  schedules.  This  allowed  stable  compressor  model  operation  with  the  nominal 
schedule  at  approximately  15°R  higher  inlet  temperatures  than  could  be  obtained  with  the 
cambered  or  axial  schedules  for  a  given  ramp  rate. 

Figure  37  shows  the  model  predicted  effect  on  compressor  stability  if  the  5,000°R/sec 
inlet  temperature  ramp  is  imposed  only  on  a  90-deg-arc  circumferential  section  of  the 
compressor  inlet.  This  case  illustrates  a  combined  effect  of  the  temperature  ramp  and  the 
associated  inlet  temperature  distortion.  Model  results  presented  in  Fig.  37  indicate  that 
compressor  surge  occurs  15.2  percent  below  the  steady-state  stability  limit  at  an  inlet 
temperature  distortion  level  of  18.4  percent. 

5.0  SUMMARY  OF  RESULTS 

Validated  TF41-A-1  HP  compressor  mathematical  models  capable  of  aerodynamic  stall 
prediction  were  developed.  The  following  summarize  the  specific  efforts  and  analyses; 


1.  A  one-dimensional,  steady-stale  model  of  the  TF41-A-1  high-pressure 
compressor  was  developed  based  on  compressor  geometry  and  stage 
characteristics  formulated  from  experimental  data  for  use  in  analyzing  steady- 
stale  compressor  phenomena  and  compressor  improvement. 

2.  The  one-dimensional,  steady-state  model  computes  compressor  surge 

pressure  ratios  within  2,4  percent  of  that  observed  experimentally  for  nominal, 
cambered,  and  axial  IGV  schedules  over  the  corrected  rotor-speed  range  of  from 
90  to  100  percent. 

3.  The  one-dimensional,  steady-state  model  computes  compressor  surge 

pressure  ratios  within  4.0  percent  of  the  experimental  stability  limits  for 
nominal,  cambered,  and  axial  IGV  schedules  with  seventh-stage  bleed-valve- 
open  operation. 
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4.  A  three-dimensional,  lime-dependent  TF41-A-1  HP  compressor  model  was 
developed  for  use  in  analyzing  transient  and  dynamic  disturbances  with  uniform 
or  distorted  inflow. 

5.  The  TF41-A-1  three-dimensional,  time-dependent  model  was  subjected  to  pure 
radial  pressure  distortion,  pure  circumferential  pressure  distortion,  and  pure 
circumferential  temperature  distortion  patterns  of  varying  severity.  Their 
influences  on  stability  were  computed.  No  specific  experimental  data  for  the 
TF41-A-1  HP  compressor  subjected  to  these  specific  distortion  patterns  were 
available.  Therefore,  comparisons  were  made  to  experimental  results  from  other 
compressors.  Trends  and  general  magnitudes  computed  by  the  model  agreed 
with  the  experimental  results. 

6.  The  TF41-A-1  three-dimensional,  time-dependent  model  was  driven  to 
instability  by  ramping  inlet  temperature  at  rates  of  2,000,  5,000,  and 
l0,000cR/sec  to  simulate  ingestion  of  hot  gases  from  gun  or  rocket  fire.  Ramp 
heights  varied  from  74  to  I08°R  with  the  2,000  and  10,000°R/sec  inlet 
temperature  ramp  rates,  respectively,  at  the  cambered  IGV  schedule. 

7.  The  aerodynamic  instabilities  that  were  encountered  during  model  loading 
simulations  occurred  in  test  surges  and  exhibited  characteristics  representing  the 
actual  physical  instability. 
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Figure  6.  Effects  of  IGV  schedule  on  TF41-A-1  HP  compressor  stability  limit  at  30,000  ft, 
Mach  number  0.7  with  seventh-stage  bleed  valve  closed. 
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Figure  7.  Effects  of  IGV  schedule  on  TF41  A  1  HP  compressor  stability  limit  at  sea  level-static 
conditions  with  seventh-stage  bleed  valve  open. 
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Figure  9.  Radial  total  pressure  profiles  with  uniform  engine  inlet  flow. 
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b.  Stages  5  through  8 
Figure  9.  Continued. 
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Figure  10.  Radial  total  temperature  profiles  with  uniform  engine  inlet  flow. 
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b.  Stages  5  through  8 
Figure  10.  Continued. 
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Figure  11.  Station  locations  {computational  planes)  for  one-dimensional,  steady-state  compressor  model. 
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Figure  12.  Elemental  control  volume  with  application  of  mass, 
momentum,  and  energy  conservation  principles. 
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Figure  13.  Typical  temperature  and  pressure  coefficients  as  functions 
of  flow  coefficient  at  constant  rotational  speed. 
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Figure  14.  Continued. 
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with  Experimental  measurements  at  30,000  ft,  MO  =  0.7. 


56 


AEDC-TR-79-39 


o» 

o 


Experimental  Data 


Sym  IGV,  deg 


NHA^,  percent 


O  18.7  ±  0.1  93.1  to  95.9 

&  11.0  ±  0.2  95.4  to  97.8 

□  10.8  ±  0.1  98.0  to  100.0 


Solid  Symbols:  Surge  Points 


b.  Cambered  IGV  schedule 
Figure  15.  Continued. 


A  EDC-TFI  -79-39 


HP  Compressor  Pressure  Ratio,  P„/P 


Experimental  Data 

Sjon  IGV,  deg  ^4'  percent 

O  4.2  +  0.1  89.9  to  91.5 

A  -0.3  ±  0.1  91 .5  to  93.3 

□  -5.5  ±  0.1  92.8  to  94.8 

k  -5.4  ±  0.1  95.1  to  96.7 


HP  Compressor  Corrected  Airflow,  w 24“^ 24^^24’  lbm//sec 


c.  Axial  IGV  schedule 
Figure  15.  Concluded. 
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Figure  16.  Model  computed  effect  of  corrected  rotor  speed  variation  for  -0.3-deg  IGV  line 
axial  IGV  schedule  at  30,000  ft,  MO  =  0.7. 
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Figure  17.  Comparison  of  model  computed  speed  characteristics  and  stability  limits 
with  experimental  measurements  with  seventh-stage  bleed  valve  open 
at  sea-level-static  conditions. 
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Figure  18.  Comparison  of  model  computed  stage  flow  coefficients 
near  compressor  stability  limit  with  experimental 
measurements,  IGV  at  full  axial  position 
<-5.0  ±  0.5  deg). 
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b.  Stages  5  through  8 
Figure  18.  Continued. 
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c.  Stages  9  through  11 
Figure  18.  Concluded. 
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Figure  19.  Comparison  of  model  computed  stage  pressure  ratios  near 
compressor  stability  limit  with  experimental  measurements, 
IGV  at  full  axial  position  (-5.0  ±  0.5  deg). 
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b.  Stages  5  through  8 
Figure  19.  Continued. 
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Figure  19.  Concluded. 
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Figure  20.  Comparison  of  model  computed  stage  temperature  ratios  near 
compressor  stability  limit  with  experimental  measurements, 
IGV  at  full  axial  position  (-5.0  ±  0.5  degt 


72 


Temperature  Ratio  Near  Compressor  Stability  Limit 


1.12 

1.08 

1.04 

1.00 

1.12 

1.08 

1.04 

1.00 

1.12 

1.08 

1.04 

1.00 

1.12 

1.08 


Alti 

Mach 

Axia! 

Sym 


Stage 


AEDC  TR-79-39 


Stage  Temperature  Ratio  Near 
Compressor  Stability  Limit 


A  ED  C-TR  -79-39 


Altitude  =  30,000  ft 
Mach  Number  =  0.7 
Axial  IGV  Schedule 


NHA£,  percent 


c.  Stages  9  through  11 
Figure  20.  Concluded. 
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Figure  23.  Schematic  of  three-dimensional,  time-dependent  solution  procedure  for  TF41-A-1  HP  compressor  model. 
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Figure  24.  Illustration  of  effect  of  incremental  averaging  on  numerical  instability. 
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Figure  25.  TF41-A-1  HP  compressor  model  loaded  to  stall,  NH/\/6"=  92.6 
percent,  IGV  =  12  deg,  no  distortion. 
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Figure  26.  TF41-A-1  HP  compressor  model  flow  breakdown  in  terms  of  computed 
stage  inlet  total  pressures  leading  to  instability,  NHA/F=  92.6 
percent,  IGV  =  12  deg,  no  distortion. 
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Figure  27.  Computed  radial  distortion  influences  on  stability, 
TF41-A-1  HP  compressor. 
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Figure  28.  Comparison  of  TF41-A-1  HP  compressor  model  radial  distortion 
computations  with  experimental  results. 
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Figure  29.  Model  computed  influence  of  circumferential  pressure 
distortion  on  stability,  TF41-A-1  HP  compressor. 


84 


Ratio  of  Stability  Limit  Pressure  Ratio  Distorted  Flow 
to  Limit  Pressure  Ratio  with  Clean  Fluw 


AEDC-TR-79-39 


Distortion  Pattern  Magnitude, 
(PH  "  PL)/Pavg’  percent 


Figure  30.  Comparison  of  TF41-A-1  HP  compressor  model  circumferential 
distortion  computations  with  experimental  results. 
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Figure  31.  Model  computed  influence  of  circumferential  temperature 
distortion  on  stability,  TF41-A-1  HP  compressor. 
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Figure  32.  Comparison  of  circumferential  temperature  distortion 
model  computations  with  experimental  results. 
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Figure  33.  Effect  of  5,000°R/sec  uniform  inlet  temperature  ramp  on  compressor  stability 
as  computed  by  TF41-A-1  HP  compressor  model. 
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TF41-A-1  HP  compressor  model  flow  breakdown  with  5,000°R/sec 
inlet  temperature  ramp,  nominal  IGV  schedule,  no  distortion. 
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Figure  36.  Effect  of  5,000° R/sec  uniform  inlet  temperature  ramp 
on  compressor  inlet  corrected  airflow  as  computed  by 
TF41-A-1  HP  compressor  model. 
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Figure  37.  Effect  of  5,000°R/sec  inlet  temperature  ramp  (imposed  on  90-deg-arc 

circumferential  section)  on  compressor  stability  as  computed  by  TF41-A-1 
compressor  model. 
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AEDC-TR-79-39 


NOMENCLATURE 

A  Area 

Az  Control  volume  area  normal  to  axial  direction 

a  Acoustic  velocity 

BVP  Bleed  valve  position 

Cd  Duct  pressure  loss  friction  coefficient 

Cp  Specific  heat  at  constant  pressure 

C„  Specific  heat  at  constant  volume 

e  Internal  energy 

F  Force  of  compressor  blading  and  cases  acting  on  fluid,  including  wall 

pressure  area  force 

f  (...)  A  function  of  ... 

g  (...)  A  function  of ... 

H  Total  enthalpy  flux 

HC  Total  enlhalpy  flux  transported  across  control  volume  circumferentially-facing 

boundary 

HP  High-pressure  rotor  ^ 

HR  T otal  enthalpy  flux  transported  across  control  volume  radially-facing  boundary 

HZ  Total  enthalpy  flux  transported  across  control  volume  axial-facing  boundary 

IGV  Inlet  guide  vane 

IMP  Impulse  function 

IMPR  Ratio  of  stage  exit  to  stage  entry  impulse  functions 

i  Axial  location  index 


93 


AEDC-TR-79-39 


j 


ip 

J 

j 

,  JP 

k 

kp 

M 

N 

NH 

NL 

P 

PD 

PR 

PS 

Q 

R 

r 

s 

T 

TR 

TS 

t 

U 


Adjacent  axial  location,  i  +  1 
Mechanical  equivalent  of  heat 
Radial  location  index 
Adjacent  radial  location,  j  +  1 
Circumferential  location  index 
Adjacent  circumferential  location,  k  +  1 
Mach  number 
Compressor  rotor  speed 
High-pressure  compressor  rotor  speed 
Low-pressure  compressor  rotor  speed 
Stagnation  (total)  pressure 
Dynamic  pressure  measurement 
Stage  total  pressure  ratio 
Static  pressure 

Rate  of  heat  addition  to  control  volume 
Gas  constant 

Radius;  coordinate  in  the  radial  direction 

Coordinate  in  the  circumferential  direction 

Stagnation  (total)  temperature 

Stage  total  temperature  ratio 

Static  temperature 

Time 

Axial  velocity 
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AEDC-TH -79-39 


v  Radial  velocity  component  in  three-dimensional  model  development 

vol  Volume 

W  Mass  flow  rate;  airflow  rate 

WB  Compressor  bleed  flow  rate 

WC  Mass  flux  across  circumferentially-facing  control  volume  boundary 

WR  Mass  flux  across  radially-facing  control  volume  boundary 

WS  Stage  shaft  work  added  to  fluid  in  control  volume 

WZ  Mass  flux  across  axially- facing  control  volume  boundary 

w  Circumferential  velocity  component  in  three-dimensional  model 

development 

X  Energy  function 

z  Axial  coordinate 

GREEK  SYMBOLS 
a  Angle  of  attack 

/?  Flow  direction  angle  relative  to  axial  direction 

7  Ratio  of  specific  heats 

A  A  difference 

5  Ratio  of  compressor  entry  total  pressure  to  standard  day,  sea-level-static 

pressure 

8  Ratio  of  compressor  inlet  total  temperature  to  standard  day,  sea-level- 
static  temperature 

\  Blade  chord-to-axial  direction  angle  or  stagger  angle 

p  Density 

<t>  Stage  flow  coefficient 
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AEDC-TR-79-39 


\J/P  Stage  pressure  coefficient 

ipT  Stage  temperature  coefficient  (stage  loading  parameter) 

SUPERSCRIPTS 
P  Pertaining  to  pressure 

T  Pertaining  to  temperature 

SUBSCRIPTS 

0,1,2  ...  Location  or  station  designators 
avg  Average  value 

exit  Pertaining  to  compressor  exit  plane 

H  High 

i  Axial  location  index 

ijk  Three-dimensional  location  index,  refers  to  value  at  location,  i,  j,  k 

ip  Adjacent  axial  location,  i  +  I 

j  Radial  location  index 

jp  Adjacent  radial  location,  j  +  1 

k  Circumferential  location  index 

kp  Adjacent  circumferential  location,  k  +  1 

L  Low 

loc  Local  value 

max  Maximum  value 

min  Minimum  value 

n  Last  location,  e.g.,  last  axial  location  in  control  volume 

p  Plus  one 
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A  ED  C-TB  -79-39 


radial  avg 
ring  avg 
wh 
z 


Average  value  along  a  radius 

Average  value  along  a  circle,  average  at  constant  radius 
Wheel 

Axial  component 


/ 


97 


